Maximizing the electric field strength in the foci of high numerical aperture optics 
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Several applications require spatial distributions of the incident electric field that maximize the electric field 
at a focal point for a given input power. The field distributions are derived for various optical systems in a direct 
way based on fundamental physical properties. The results may prove useful for a wide range of applications, 
e.g., microscopy, scattering experiments or excitation of single atoms. For commonly used distributions - fun- 
damental Gaussian modes and doughnut modes - we give the upper bounds of the achievable field amplitudes. 



I. INTRODUCTION 

Recently, there has been increased interest in the interac- 
tion of single matter systems like atoms, molecules or ions 
with weak light beams or single photons in free space. The 
topical spread covers scattering by single two-level systems 
(HI |2] |3] m [51 as well as the absorption of single photons by 
single atoms f6ll7]f8^|. Also other applications relying on high 
numerical aperture (NA) optics have raised increasing atten- 
tion, involving parabolic mirrors El[l3 as well as lenses 1 11 1. 

It is common to all these applications that they require the 
smallest possible focal volume or - in other words - that the 
power incident onto the focusing optics is transformed in such 
a way that the electric energy density in the focus is maxi- 
mized. Thus, one wants to maximize the amplitude of the 
electric field in the focus. Although it is well known since 
quite some time that the optimum performance is achieved 
by electric dipole radiation lfT2l [T3l [141 and also some newer 
publications hint in this direction fT5"T6l, this recipe is fol- 
lowed only in few papers in a consequent manner dS] [6) [TJ . 
There are also publications in which some sort of 'optimiza- 
tion procedure' is followed: In Ref. f4] the beam width of a 
circularly polarized fundamental Gaussian is varied in order 
to maximize the strength of the circularly polarized field com- 
ponent in the focus of an aspheric lens. The authors of Ref. 
|[TOl follow a similar procedure in varying the beam waist of a 
radially polarized doughnut beam during minimization of the 
focal volume in the case of a parabolic mirror In both cases, 
the optimization procedure is equivalent to maximizing the 
overlap of the incident radiation with dipole radiation. In the 
formercase |4| the Gaussian fundamental mode is matched to- 
wards the distribution that creates the field of a circular dipole 
after refraction by the focusing lens. In the latter case fTU\ the 
doughnut mode is adapted such that it best overlaps with the 
field of a linear dipole that has its axis parallel to the optical 
axis of the parabolic mirror Finally, in a recent paper the field 
distribution that maximizes the focal field component paral- 
lel to the optical axis of a high NA lens is directly calculated 
by a variational method L 11 J . However, with the exemption 
of Ref. ifTOll no direct and explicit relation between the found 



optimal radiation modes at the input and a dipole radiation 
pattern is given in these publications. 

Based on earlier results obtained by others |[T2l[T3l we fol- 
low a different approach in what follows here. Knowing that 
it is only electric dipole radiation that contributes to the focal 
electric field (see, e.g., Refs. |[T2l[r6!l ), we calculate the distri- 
bution of the field incident onto the focusing device that gives 
the desired dipole radiation pattern after reflection/refraction 
by the focusing optics, as it was done for a parabolic mirror 
and linear dipole radiation in Ref. |6 1. This approach is on the 
one hand the most direct one. On the other hand, it is concep- 
tually simple in comparison to other methods. Besides tem- 
poral/spectral issues not treated here, this ansatz is the same 
as the time reversal argument given in Refs. ||6l|7l (see also, 
e.g., Ref. flTl for time reversal focusing of microwaves and 
Ref |[T8l for acoustic waves). 

Some of the results obtained here can in parts also be found 
in other publications |[T3l [T9l [20l [2T 1. However, a compre- 
hensive study that covers all practically relevant cases is still 
lacking. Therefore, this paper is intended to deliver a recipe 
for maximizing the electric field in focusing applications, in 
particular the excitation of single atoms as envisaged in Refs. 
iH [4] [5] [71 • In the next section, we will emphasize the impor- 
tance of a large solid angle of illumination. In Sees. Ill and 



[IV[the optimal intensity and polarization distributions will be 
derived for the cases of illumination by means of a parabolic 
mirror, a lens fulfilling the sine condition and an ideal thin 
lens. Since the experimental realization of these distributions 
may turn out to be rather involved, the overlap of 'standard' 
radiation modes with the ideal ones is treated in Sec. jV] 

Unlike in many other publications where the field in the fo- 
cal region is treated, the discussion will be restricted to the 
electric field in the very focal point itself. For the main appli- 
cations in mind here - scattering by single atoms or absorption 
of single photons - this constitutes a well justified restriction, 
since it is only the field at the location of the atom that matters 
(see also Ref |4 , TSl). Furthermore, this restriction leads 
to several simple analytical results. 



n. INFLUENCE OF THE SOLID ANGLE 
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Before starting with the calculations, we make some notes 
on the terminology used in this paper. For the quantum-optical 
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applications mentioned above, the field parallel to an atomic 
dipole has to be maximized. In these cases, the 'dipole axis' 
is given by the quantization axis of the scenario at hand. In a 
classical sense, the term 'dipole axis' designates the direction 
from which the polar angle d is measured that parameterizes 
the angular dipole radiation patterns (see also Fig. [T]l. How- 
ever, one may also think of applications in which no particu- 
lar emitter or source is present at the focus but one wants to 
maximize the electric field vector that points into a certain di- 
rection. This field vector can be thought of as being parallel to 
some imagined dipole moment. Thus, we will adopt the ter- 
minology of Ref |5 1 and often speak of 'virtual dipoles', i.e., 
dipoles that would produce a certain kind of radiation pattern 
but are not present in reality. 

The term 'solid angle' usually designates the integral over 
the covered polar angle and the covered azimuthal angle. 
Here, we use this term in a slightly different context. It des- 
ignates the covered solid angle weighted by a certain dipole 
emission pattern: Ail = j D{d) sin ■&<!■& dip, W\t\\ DTj{-d) = 
sin^ t9 for a Unear dipole or D„j_ (i?) = (1 + cos^ /2 for a 
circular dipole f22], respectively. Thus, the maximum achiev- 
able, weighted solid angle is Afimax — Stt/3 instead of in in 
the usual sense. 

We now derive an analytic expression that highlights the in- 
fluence of the covered solid angle. The electro-magnetic field 
can be decomposed into its different multipole components. 
This has been exploited in Ref |20| to calculate the field in 
the focal region for a given input field. However, since we are 
only interested in the field at the focus, which we take to co- 
incide with the origin of our coordinate system, we only need 
to take care of the dipole components. Furthermore, it has al- 
ready been pointed out in Ref. 1 19 1 that each kind of dipole 
radiation creates predominantly its respective state of polar- 
ization at the focus/origin. This is exactly true for dipoles 
oriented along the optical axis of the focusing system. There- 
fore, it is obvious that one only needs to create the radiation 
pattern of a single kind of dipole for most practical cases. 

In order to calculate the (time independent) field in the fo- 
cus, we use the Debye-integral in the form given in Ref. 
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where E/(i?, p) is the field on the focal sphere for a focal 
length /, ip is the azimuthal angle and A the wavelength of the 
light to be focused. Furthermore, a term exp(ikr) has already 
been dropped since we treat only the case r = 0. 

If only a single kind of dipole radiation is incident onto 
the focus, the field on the focal sphere can be written as 
E/ = ipdPd, where i/^d is a proportionality factor and is 
the polarization vector of dipole radiation in the far field with 
d — Tr,a± for a linear or circular dipole, respectively. The 
polarization vectors in spherical coordinates are given by [24J 



Pti — — sin 'de^ 
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with 6,5 being the polar and azimuthal unit vectors. 

Next, we require that independent of the covered solid an- 
gle the power P incident onto the focus is a constant (this 
would correspond to the case of, e.g., a single photon that is 
focussed with lenses of different numerical aperture but al- 
ways covers the whole lens cross section). The total power 
radiated through the surface of the focal sphere can be written 
as 
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where the integration limits are given by the intervals of the 
azimuthal angle and the polar angle covered by the incident 
light. The above expression is solved for ipd'- 
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where Afi^ is the solid angle weighted by the dipole radiation 
pattern: 



Afi^ 

for a linear dipole and 
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for a circular dipole. 

Inserting the corresponding polarization vector and the 
factor Ipd into Eq. [Tjand performing the integration for a linear 
dipole as well as for circular dipoles, one obtains a general 
relation for the field component parallel to the virtual dipole: 



E<i(0) 



(8) 



and 



This emphasizes that in order to maximize the focal field com- 
ponent parallel to a given virtual dipole one has to (i) put as 
much of the incident power as possible into the corresponding 
dipole mode and (ii) maximize the solid angle covered by the 
incident radiation. It should be pointed out that the resulting 
field strength does not depend on the dipole type. 

We note that for full solid angle illumination the field given 
by Eq. [8] delivers exactly the maximum possible electric en- 
ergy density given by Bassett lfT2l . Furthermore, E q. [8| consti- 
tutes the generalization of Eqs. 8 and 10 of Ref. fill], where 
the case of a linear dipole oriented along the optical axis of an 
objective obeying the sine condition is treated. However, the 
connection to dipole radiation was not established there. 



III. OPTIMUM IRRADIANCE DISTRIBUTIONS 



Pct± — — ^ {cos 11) [cos (fi ±i sin p]e^ + [±i cos If — sin ip]e^) 
v2 

(3) 



Next, we want to derive the ideal radiation patterns incident 
onto the focusing device which are transformed into a dipole 
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FIG. 1 : Relation of the angles a and i9 to the Cartesian coordinates: 

(a) dipole/quantization axis (gray arrow) parallel to the optical axis; 

(b) dipole/quantization axis perpendicular to the optical axis. 



wave that moves towards the focus. For all devices treated 
here, we assume that they operate perfectly, i.e., the wave in- 
cident onto the focus is assumed to have a uniform phase af- 
ter reflection/refraction off the focusing device (see also Ref. 
ifTOl ). For example, a parabolic mirror might exhibit (small) 
deviations from the perfect parabolic shape. Furthermore, the 
phase shift upon reflection is different for different angles of 
incidence. We presuppose that such effects are compensated 
for by, e.g., appropriate correcting elements. The production 
of such elements and the measurement of imperfections might 
turn out to be a challenging but feasible experimental task 

The method is simply to start with a dipole wave which vir- 
tually emerges from the focus and trace it back through the op- 
tical element under consideration. We have chosen to formu- 
late the ideal radiation patterns in terms of irradiances instead 
of field vectors, since the former quantity is the one usually 
measured in an experiment. Furthermore, in the derivation of 
the apodization factors caiTied out below, energy conservation 
considerations require the use of irradiances and radiant inten- 
sities anyway. The ideal vectors fields can be constructed from 
the irradiance patterns and the polarization patterns derived in 
the next section. 

The derivation of the intensity patterns is pursued in the 
same way as it was done in Ref. f6l for a parabolic mirror and 
a linear dipole oriented parallel to the optical axis of the mir- 
ror. We designate by a the angle that is enclosed by the optical 
axis of the focusing element and the propagation direction of 
a single ray which hits the optical element at a distance r from 
the optical axis. Here, we chose the optical axis to coincide 
with the z-axis (cf. Fig. [T]l. Energy conservation demands 
that t6J 



T(a) sin ado. = I{r)r dr 



which delivers 



with A{r, a) 



I{r) ^ T{a)-A{r,a) 

sin a(r) da{r) 
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being an apodization factor 
T{a) is the radiant intensity, i.e., the light power emitted into 
an infinitesimal solid angle. I{r) is the irradiance of the in- 
cident plane wave, i.e., the light power incident onto an in- 
finitesimal surface area. 

In what follows, we express the apodization factor as a 
function of the Cartesian coordinates x, y in the plane per- 



pendicular to the optical axis: A{x,y) = A{r{x,y)) with 
r — \J x^ + y^- In the case of a parabolic mirror, a is given 
through the relation [6\ 
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which with the help of some algebra leads to 
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For an aplanatic lens that fulfills the sine condition, i.e., a lens 
where the refracted rays of same phase seem to emerge from 
a spherical surface around the focus, one has 



r 
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resulting in 
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(13) 
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As a third focusing element, we consider an ideal thin lens, 
i.e., a lens for which the refracted rays seem to emerge from a 
plane perpendicular to the optical axis. A practical realization 
of such a thin lens would be an idealized diffractive optical 
lens. One has 



tan ~ J 



which delivers 



/^(fl + f^ + l) 
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(15) 



(16) 



Next, for every point {x,y) in the plane perpendicular to 
the optical axis we determine The radiant intensity along 
this direction is given by the angular dipole radiation pattern 



(17) 



where Dq is a proportionality constant. Finally, T(i9) is also 
expressed in terms of x,y and multiplied with the apodization 
factor A{x, y) yielding the irradiance I{x, y). 

If the dipole is oriented parallel to the optical axis, a co- 
incides with -d (see Fig. [T^). Hence, the ideal irradiation is 
rotationally symmetric with respect to the optical axis and the 
angular dipole emission patterns given above can be written 
as 



'T,,a±{a). Using Eqs. 



11 



13 



and 



15 



^ 

spectively, and applying several trigonometric relations one 

obtains the ideal irradiance distributions. They are summa- 
rized in Tab. [l] and plotted in Fig. [2] We have also intro- 
duced normalized coordinates X ~ X / f ,Y — y/f, R = r/ f. 
All remaining dimensional quantities have been lumped into 
a proportionality constant /q. 
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TABLE I: Optimum irradiance distributions I{R) for the creation of radiation patterns of dipoles oriented parallel to the optical axis of the 
focusing device. R = r/ f is the radial coordinate transverse to the optical axis given in units of the focal length. The proportionality constant 
Jo contains all dimensional quantities. 
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FIG. 2: Ideal irradiance distributions for dipoles oriented parallel 
to the optical axis of the focusing device: (a) parabolic mirror, (b) 
aplanatic lens fulfilling the sine condition, (c) ideal thin lens. Solid 
(dashed) lines denote the irradiance for a linear (circular) dipole ra- 
diation pattern. 



FIG. 3: Ideal irradiance distributions for dipoles oriented orthogonal 
to the optical axis of the focusing device: (a) parabolic mirror, (b) 
lens fulfilling the sine condition, (c) ideal thin lens. Left panel: linear 
dipole; right panel: circular dipole. The gray levels designate values 
oil I la ranging from (black) to 1 (white) in (a,c) and from (black) 
to 3 (white) in (b). 



The formula obtained here for the irradiance in the case 
of an aplanatic lens fulfilling the sine condition and a linear 
dipole corresponds exactly to the square of the formula for 
the electric field amplitude in the entrance pupil of the lens as 
it was found by means of an optimization method in Ref. ifTTI . 

Next, we turn to the case of the virtual dipole being ori- 
ented perpendicular to the optical axis of the focusing ele- 



ment. Without loss of generality we choose the virtual dipole 
axis (the direction i9=0) to coincide with the y-axis (see Fig. 
[TJi). An immediate consequence of the dipole axis being or- 
thogonal to the optical axis is that the ideal intensity distribu- 
tions have no longer rotational symmetry with respect to the 
optical axis. Furthermore, the angles a and -d do not coincide 



any more. Thus, the relations given by Eqs. 11 13 and 15 can 
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not be used to express the angular emission patterns D-^ ,j± {"&) 
in terms of the transverse coordinates. Instead one has to use 
in the case of the parabolic mirror the relations 



y 



,z = f- 
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(18) 



This is valid if the focus of the parabolic mirror is located in 
the origin of the coordinate system. For a lens obeying the 
sine condition the relation between -d and x, y reads 



cos d = - 
] 

and for the ideal thin lens one has 



tan-i? 



y 



(19) 



(20) 



Again, the focus is located in the origin of the coordinate sys- 
tem. Using some trigonometric relations and multiplication 
with the corresponding apodization factors leads to the ideal 
irradiance distributions given in Tab. [lljand displayed in Fig. 



IV. OPTIMUM STATES OF POLARIZATION 



Again, we begin the discussion for a virtual dipole parallel 
to the optical axis. For a linear dipole, the polarization pattern 
can be obtained immediately from the fact that the polariza- 
tion vector of a linear dipole is oriented along in the far 
field. This results in the well known radial polarization pat- 
tern ll6l [T4l[T5]| and, of course, coincides with the result found 
recently by optimization ifTTI . The polarization pattern is the 
same for all of the devices discussed here. 

The situation is slightly more complex in the case of a cir- 
cular dipole. The polarization vector for a certain is given 
by Pctj_ in Eq. [31 Using the relations given in Ref. p26l| the 
normalized Stokes parameter S3 / So describing the degree of 
circular polarization is readily determined to be 



S'scr, ±2 cost? 



1 + cos^ d 



(21) 



where Sq corresponds to the total intensity and S3 to the inten- 
sity that would be measured after passing the beam through a 
circular polarizer With S3 /So = sin(2x) giving the so called 
ellipticity angle |27|, the ellipticity as a function of {) is 
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(22) 



where x<j± = ±vr/4 corresponds to circular states of polariza- 
tion. This shows that the ideal polarization pattern for circular 
dipoles exhibits a varying ellipticity as a function of the radial 
distance to the optical axis. 
Using Eqs. 



11 13 and 15 to express i3 through r, one ob- 
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FIG. 4: Ellipticity of the ideal states of polarization for a circular 
dipole ((T+ , a sign change gives the values for a_ ) oriented parallel to 
the optical axis of the focusing device. Solid line: parabolic mirror; 
dashed line: aplanatic lens; dash-dotted line: ideal thin lens. 



17 7 ^'^rlS 
^ -1 7 r 

7 ^ -J f J' ^ ^ 

A -i 

I /I ^ //^ ~ 
jfi /f // ✓ — J 
« / / -7 -7 *!. 

r r 1 h 





-7 1 ir r 

^ /I / 
c-,' 4 ^ A 

J' TT '1 y\ 7 

' .7 r -7 



tains the ideal state of polarization for the different focusing 



FIG. 5: State of polarization of the ideal radiation pattern for the 
case of a parabolic mirror and a circular dipole ((7+) with the dipole 
axis oriented parallel to the optical axis of the mirror. The vectors in 
gray color denote the linearly polarized components with 7r/2 phase 
shift. The dashed line denotes the boundary NA=1. 



devices. The results are displayed in Fig. |4] Note that in 
the case of the parabolic mirror there is a sign change of the 
ellipticity at the radius r — 2/, which corresponds to a half 
opening angle of the parabolic miiTor of 90° . This result ex- 
plains the finding of Ref. |10|(Fig. 2) that an increase of the 
half focusing angle beyond 90° does not lead to a substantial 
decrease of focus size for an incident wave that has a circu- 
lar state of polarization of same helicity over the whole beam 
cross section. 

The spatial polarization pattern in the entrance plane of 
the focusing devices is calculated exemplarily for a parabolic 
miiTor First, p^^ is transformed into a unit vector \>a^ by 
multiplication with \/2/^/l + cos^ d. Then, the polarization 
pattern after reflection off the parabolic miiTor is obtained as 
follows: The state of polarization p^ ^^ of the virtual dipole, 
formulated in terms of the spherical unit vectors ef)_^, is trans- 
formed into the Cartesian coordinate system: p^.i^ Px,y.z- 
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TABLE II: Optimum irradiance distributions I{X, Y) for the creation of radiation patterns of dipoles oriented orthogonal to the optical axis 
of the focusing device. X = x/f and Y = y/ f are the coordinates transverse to the optical axis given in units of the focal length. The 
proportionality constant Jo contains all dimensional quantities. 
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Then, the ray for each emission direction is traced from the 
virtual dipole towards the surface of the parabolic mirror The 
polarization vector of the ray after reflection at the mirror sur- 
face is calculated via the following equation: 



P = 2(p^,j,,^ • N)N - p. 



x,y,z 



(23) 



with the surface normal of the mirror given by 
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(24) 



Again, the focus is chosen to be at the origin of the coordinate 
system. For a± light, this procedure results in 
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(25) 



This pattern is depicted in Fig. [5] Up to the boundary given 
by \/X^~+Y^ = 2 (or numerical aperture NA=1), the polar- 
ization pattern is qualitatively similar for the other focusing 
devices treated here. However, one has to keep in mind that 
the relations between -d, ip and the Cartesian coordinates in the 
entrance plane of the devices are different. Thus, the pattern 
displayed inside the circle marking NA=1 in Fig.|5]will be dis- 
torted or stretched/compressed to some amount, respectively. 

If the dipole axis is orthogonal to the optical axis, the pro- 
cedure is the same except that in this case the direction -d — 
coincides with the positive y-axis. This is handled by inter- 
changing the second with the third component of the unit vec- 
tors of the spherical coordinate system. Using Eqs. |2]and 
[3] one obtains the polarization patterns 



v/(x2 + r2 + 4)2 



i6y2 




for a linear dipole and 
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(27) 



for the circular dipoles. The two polarization patterns are il- 
lustrated in Fig.|6] 



V. OVERLAP OF PRACTICAL DISTRIBUTIONS WITH 
THE IDEAL ONES 

Some of the radiation patterns derived above - namely 
in the case of the dipole axis parallel to the optical axis - 
exhibit obvious similarities with radiation modes which are 
commonly used in experiments. In the case of a parabolic 
mirror and the ideal thin lens, the irradiance distributions for 
a linear virtual dipole look close to the well known doughnut 
modes. In the case of circular virtual dipoles there is a strong 
similarity with fundamental Gaussian distributions. 

In order to estimate and quantify the strength of the simi- 
larities, we will calculate the overlap of these modes with the 
ideal radiation distributions. The discussion is limited to the 
case that the dipoles are oriented parallel to the optical axis. 
We begin with defining the overlap on the focal sphere as 



(28) 



77 = 27r / E* (i^) • E/ sin i9 M 



where E^ is the ideal field and Ej the incident field. The 
integration is performed over the polar angle covered by the 
focusing optics. For the sake of simplicity, we make the ap- 
proximation that the polarization vector of the incident field is 
everywhere parallel to the polarization vector of the ideal dis- 
tribution. We can then write the overlap in terms of the field 
moduli or the radiant intensities T{'d), respectively, as 



11 = 
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Identifying -d with a in Eq. |9]and taking the square root of Eq. 
|9]we have 
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FIG. 6: State of polarization of the ideal radiation patterns for the 
case of a parabolic mirror and the dipole axis oriented orthogonal to 
the optical axis of the mirror: (a) linear dipole, (b) circular dipole 
In (b) the vectors in gray color denote the linearly polar- 
ized components with n/2 phase shift. The dashed line denotes the 
boundary NA=1. 



2tt J p*^-EfSmi}di}^i;d-And , (33) 

since tlie different kinds of dipole fields are orthogonal for 
full azimuthal angle integration. Dividing both sides by 
(27r / |E/|2 sin i9 di? • 27r / |pd|2 sin I? dd)^/^ gives 



2n /lE/Psinf^dt? 



(34) 



with riN,d the normalized overlap on the focal sphere calcu- 
lated for a specific ideal dipole radiation of kind d. For any 
given input power, the field in the focus parallel to a certain 
virtual dipole is maximized when the complete incident power 
is contained in the corresponding dipole wave, i.e., t/j^ is max- 
imized: Ef = V-'d.maxPd- This mcans that 
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71 sin Z9 dt? = Vd.max ' 
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Inserting this into Eq. [34]finally reveals that 
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lpd.,1 



(35) 



(36) 



Thus, the fraction of the achievable field modulus in the focus 
with relation to the maximum possible one is directly propor- 
tional to the overlap rjN,d on the focal sphere or the overlap 
in the entrance plane of the focusing optics, respectively. In 
other words, if the incident wave is not a pure dipole wave, 
the right hand side of Eq. [8] has to be multiplied by a factor 

VN^d- 
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d ■ VN,d 



(37) 



This is the (expected) result that the overlap on the focal 
sphere is the same as the overlap in the entrance plane of 
the focusing optics. Finally, we normalize the overlap to the 
power contained in the ideal and the incident field and arrive 
at 



Vn = 



^Jh{r)rdr-^pf{r), 



(31) 



dr 



Next, the normalized overlap is related to the field ampli- 
tude in the focus. We expand the field on the focal sphere in 
terms of the electric dipole components as 



E 



7 



(32) 



Magnetic dipole terms and higher order multipole terms are 
neglected since they do not contribute to the electric field in 
the focus anyway ||T2I . Taking the integral J dipdd sin iJpJJ on 
both sides, where d= ■K,a± and is given by Eqs. |2]and[3] 



Figure |7] shows some examples for illustration. In the case 
of a parabolic mirror with a half opening angle of 1 34° as it 
is used in the experiment described in Refs. ||6] [71, a radi- 
ally polarized doughnut mode i^yi{r) ^ r ■ e"*" /™ , w is 
the beam radius) best matches the field distribution for a lin- 
ear dipole at a beam radius of w/ f —2.26. There, the over- 
lap is riN^TT =0.982. The half opening angle of 134° corre- 
sponds to Afi^ =0.94 x 87r/3. Thus, the field amplitude par- 
allel to the virtual linear dipole would be 0.95 times the one 
obtained with the ideal radiation pattern and a parabolic mir- 
ror of infinite dimensions. This has been confirmed in simula- 
tions. In the case of a circular dipole, the maximum overlap of 
VN,tT± =0.978 is achieved for a Gaussian fundamental mode 
(V^r(r) - e-'''/"'') with abeam waist of w// =3.14. With 
AOct± =0.80x Stt/S one achieves a field amplitude parallel to 
the virtual circular dipole of 0.87 times the maximum possible 
one. 

For an ideal thin lens (Fig. [T]^) with NA=0.95 the opti- 
mum overlap of i]n,tt =0.983 is reached at w// =1.84 (linear 
dipole) and the maximum value of 7]N,a± =0.985 reached at 
w/f =1.94 (circular dipole). With Afl^ =0.27x87r/3 and 
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0! 




1 2 3 4 1 2 3 

w/f w/f 



FIG. 7: Normalized overlap of the ideal field distributions with com- 
mon radiation modes as a function of the modes beam radius w in 
units of the focal length /. The virtual dipoles are oriented parallel 
to the optical axis, (a) Parabolic mirror with a half opening angle of 
134°. (b) Ideal thin lens with a numerical aperture of 0.95. Solid 
lines: Overlap of the ideal field for linear dipoles with a radially po- 
larized doughnut mode. Dashed lines: Overlap of the ideal field for 
circular dipoles with a fundamental Gaussian mode that has an ideal 
state of polarization. 




0.25 0.5 0.75 1 

half opening angle (rad/n) 



FIG. 8: Solid angle Af2^/ (Stt/S) (dash-dotted line), maximum pos- 
sible overlap riN,7r with a linear dipole wave (dashed line) and field 
amplitude in the focus in units of the maximum possible one 
(solid line) as the function of the half opening angle of a parabolic 
mirror. For further explanations see text. 



TABLE III: Maximum electric field amplitudes achievable by use 
of a radially polarized doughnut mode and a fundamental Gaussian 
mode with the state of polarization matching a circular dipole. The 
values are given relative to the ones obtainable with the ideal irra- 
diance distributions of Fig. |2|a,c). The virtual dipoles are oriented 
parallel to the optical axis. 





doughnut mode 


fundamental Gaussian 


parabolic mirror 


0.958 


0.919 


ideal thin lens 


0.575 


0.628 



field in the focus parallel to a specific virtual dipole is max- 
imized by maximizing the solid angle, as long as the simple 
mode is polarized parallel to the ideal polarization across the 
whole beam cross section. This is illustrated in Fig. [8] for a 
linear dipole and a parabolic mirror (in the case of a circu- 
lar dipole and/or a thin lens the results are qualitatively simi- 
lar). For each (half) opening angle of the mirror the maximum 
overlap is obtained from a curve like the one in Fig. |7ja). The 
resulting field is then obtained from Eq. [37] by inserting the 
overlap value and the corresponding square root of the cov- 
ered solid angle into this equation. 

From the asymptotic behavior of the plot for the field am- 
plitude in Fig. [8]one can extract the maximum field amplitude 
one could achieve with a doughnut mode of suitable polariza- 
tion by focusing the mode with a parabolic mirror Likewise, 
one obtains the same values for circular dipoles (i.e., overlap- 
ping with a Gaussian mode of proper polarization) and/or an 
ideal thin lens. The values are listed in Tab. [Ill] 

One could calculate the overlaps also in the case of a lens 
obeying the sine condition. However, due to the divergence 
of the ideal irradiance values towards the boundary of the lens 
aperture the calculation of the overlap with the simple modes 
discussed above does not appear to be meaningful. Similar 
arguments apply when the virtual dipole is oriented orthogo- 
nal to the optical axis. Since the ideal irradiance patterns have 
no rotational symmetry, the overlaps will always be smaller 
than their counterparts for the dipole axis being parallel to the 
optical axis. 



Ario-j- =0.38x87r/3 one achieves 0.51 and 0.61 times the 
field amplitudes of the maximum possible ones, respectively. 

Please note that for the parabolic mirror as well as for the 
ideal thin lens the optimum beam waist of the doughnut is 
such that the maximum of the irradiance distribution is located 
at a somewhat larger radial position than the maximum of the 
radiation patterns of a linear dipole that are displayed in Fig. 
|2|a,c). Furthermore, the values for the circular dipoles given 
above have to be considered as upper limits in the sense that 
the ideal state of polarization for a circular dipole is difficult 
to obtain in practice. 

If the solid angle covered by the focusing optics is in- 
creased, the maximum achievable overlap of the 'simple' 
modes with the ideal distributions decreases. However, this 
is overcompensated for by an increase in solid angle, i.e., the 



VI. CONCLUDING REMARKS 

We have derived the ideal radiation patterns with which one 
achieves the maximum possible field strengths for several fo- 
cusing devices. Besides using these radiation patterns, it is of 
utmost importance to maximize the covered solid angle, as it 
is shown by Eqs. [8]and 37 Thus, the maximum possible elec- 
tric field strength is achieved by covering the full solid angle. 
Also, the maximum possible effect in the interaction of light 
with single atoms will also only be found for full solid angle 
coverage, since scattering as well as absorption of photons is 
proportional to the square of the electric field amplitude at the 
location of the atom |26|. This finding is in contrast to the 
statements made in other publications, where it was claimed 
that the power scattered by a single atom is already maximized 
for half solid angle illumination |4, 5 1. According with the pa- 
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pers of Van Enk and Kimble lfT6l l28l . we want to state once one has to shape the radiation incident onto the atom to re- 
more that for optimum interaction of Hght with single atoms semble an electric dipole wave. 
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